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Early  prediction  of  the  adverse  outcomes  associated  with  heat  stress  is 
critical  for  effective  management  and  mitigation  of  injury,  which  may 
sometimes  lead  to  extreme  undesirable  clinical  conditions,  such  as 
multiorgan  dysfunction  syndrome  and  death.  Here,  we  developed  a 
computational  model  to  predict  the  spatiotemporal  temperature  distri¬ 
bution  in  a  rat  exposed  to  heat  stress  in  an  attempt  to  understand  the 
correlation  between  heat  load  and  differential  organ  dysfunction.  The 
model  includes  a  three-dimensional  representation  of  the  rat  anatomy 
obtained  from  medical  imaging  and  incorporates  the  key  mechanisms 
of  heat  transfer  during  thermoregulation.  We  formulated  a  novel 
approach  to  estimate  blood  temperature  by  accounting  for  blood 
mixing  from  the  different  organs  and  to  estimate  the  effects  of  the 
circadian  rhythm  in  body  temperature  by  considering  day-night  vari¬ 
ations  in  metabolic  heat  generation  and  blood  perfusion.  We  validated 
the  model  using  in  vivo  core  temperature  measurements  in  control  and 
heat-stressed  rats  and  other  published  experimental  data.  The  model 
predictions  were  within  1  SD  of  the  measured  data.  The  liver  dem¬ 
onstrated  the  greatest  susceptibility  to  heat  stress,  with  the  maximum 
temperature  reaching  2°C  higher  than  the  measured  core  temperature 
and  95%  of  its  volume  exceeding  the  targeted  experimental  core 
temperature.  Other  organs  also  attained  temperatures  greater  than  the 
core  temperature,  illustrating  the  need  to  monitor  multiple  organs 
during  heat  stress.  The  model  facilitates  the  identification  of  organ- 
specific  risks  during  heat  stress  and  has  the  potential  to  aid  in  the 
development  of  improved  clinical  strategies  for  thermal-injury  pre¬ 
vention  and  management. 

computational  modeling;  core  temperature;  finite  element  method; 
multiorgan  dysfunction  syndrome;  radiotelemetry 


heat  stress  is  characterized  by  an  increase  in  body  tempera¬ 
tures  from  exposure  to  hot  weather  conditions  (6).  Severe  heat 
stress  leads  to  heat  stroke,  which  may  trigger  a  complex  cascade 
of  events  emanating  from  thermoregulatory  failure  to  systemic 
inflammatory  response  and  eventually  to  multiple  organ  dysfunc¬ 
tion  syndrome  (MODS)  (6,  24,  55).  It  is  a  life-threatening  condi¬ 
tion  and  may  also  lead  to  increased  risk  of  mortality  in  the  long 
term  (3,  79).  Military  personnel  (4)  as  well  as  athletes  (11) 
undergoing  strenuous  physical  activity  in  hot  and  humid  environ- 


Address  for  reprint  requests  and  other  correspondence:  J.  Reifman,  DoD 
Biotechnology  High  Performance  Computing  Software  Applications  Institute, 
Telemedicine  and  Advanced  Technology  Research  Center,  U.S.  Army  Medical 
Research  and  Materiel  Command,  MCMR-TT,  504  Scott  St.,  Fort  Detrick,  MD 
21702-5012  (e-mail:  jaques.reifman.civ@mail.mil). 


ments  are  at  great  risk  for  exertional  heat  stress,  which  may  occur 
in  healthy  young  adults  (61,  89).  For  example,  in  the  United  States 
military,  14,577  cases  of  heat-related  injuries  were  reported  be¬ 
tween  2007  and  2011,  which  resulted  not  only  in  significant 
lost-duty  days  and  associated  medical  costs  but  also  in  life- 
threatening  conditions  (4,  10).  One  potential  strategy  to  reduce  the 
high  incidence  of  heat-stress  cases  and  minimize  injury  progres¬ 
sion  is  to  establish  the  ability  to  accurately  predict  impending  rises 
in  organ  temperatures  (37,  38).  In  fact,  the  need  for  the  early 
prediction  of  impending  temperature  rises  and  the  adverse  effects 
associated  with  heat  stress  has  been  recognized  as  an  integral  part 
of  effective  injury  management  and  mitigation  strategies  (6, 
39,  55). 

Injury  progression  from  the  initial  thermal  insult  to  MODS 
and  potential  mortality  is  a  complicated  process,  and  in  vivo 
animal  studies  have  been  frequently  used  to  elucidate  the 
associated  pathogenesis  and  pathophysiology  (7,  16,  54,  88).  In 
such  studies,  heat  stress  is  generally  induced  by  placing  the 
animal  in  a  temperature-controlled  chamber  where  its  core 
temperature  is  remotely  monitored  using  telemetry  probes 
surgically  implanted  into  the  abdominal  cavity  or  via  measure¬ 
ments  of  rectal  or  skin  temperature  using  thermocouples  (7,  16, 
54,  88).  However,  in  vivo  monitoring  in  these  studies  is  limited 
to  transient  temperature  measurements  at  a  single  location  in 
the  animal,  failing  to  provide  information  about  the  tempera¬ 
ture  distribution  and  presence  of  “hot  spots”  in  the  various 
organs.  The  evaluation  of  temperature  distribution  and  identi¬ 
fication  of  hot  spots  are  important  as  the  order  and  severity  of 
organ  dysfunction  in  MODS  is  dependent  on  the  inflammatory 
response,  which  in  turn  depends,  at  least  partially,  on  the 
thermal  injury  (6,  24,  80).  Because  differences  in  anatomy  and 
heat  transfer  properties  can  lead  to  distinct  thermal  injury  in 
different  organs,  invoking  differential  tissue  damage  and  organ 
dysfunction  (42,  53,  80),  it  is  critical  to  spatially  monitor  the 
temperature  in  various  organs  to  obtain  greater  insights  on 
disease  progression.  In  the  absence  of  experimental  methods,  a 
computational  model  based  on  a  realistic  description  of  an 
animal’s  anatomy  and  capable  of  emulating  the  underlying 
thermoregulatory  mechanisms  would  be  useful  for  quantifying 
the  spatial  distribution  and  temporal  dynamics  of  temperatures 
in  the  entire  body.  Such  a  model  would  facilitate  the  prediction 
of  the  time  and  extent  of  organ  damage  as  a  function  of  thermal 
stress. 

Even  under  normothermic  conditions,  visceral  organs  have 
different  temperatures,  which  depend  on  the  heat  produced  by 
cell  metabolism  and  heat  exchange  due  to  blood  perfusion  in 
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the  particular  organ  (15).  Figure  1  shows  the  various  mecha¬ 
nisms  of  heat  transfer  inside  the  body.  Conduction  and  con¬ 
vection  due  to  blood  perfusion  facilitate  heat  transfer  between 
the  various  organs.  Heat  is  lost  by  convection,  radiation,  and 
evaporation  of  sweat  or  saliva  (in  animals  that  do  not  have 
sweat  glands)  from  the  skin  surface  at  normothermia.  Conse¬ 
quently,  owing  to  these  heat-loss  mechanisms,  the  skin  is 
typically  at  a  lower  temperature  compared  with  the  core  (15, 
33).  During  heat  stress,  external  environmental  conditions 
and/or  exertion  lead  to  an  increase  in  body  temperature,  trig¬ 
gering  a  complex  physiological  response  that  aims  to  restore 
normothermia  (24).  An  increase  in  body  temperature  during 
heat  stress  induces  an  increase  in  cardiac  output,  skin  blood 
perfusion,  and  perspiration/salivation  to  facilitate  heat  loss  (6). 
The  blood  vessels  in  the  visceral  organs  vasoconstrict,  whereas 
the  skin  blood  vessels  vasodilate  to  accommodate  blood  redis¬ 
tributed  from  the  core  to  the  body  surface  by  the  cardiovascular 


Model  Formulation: 

Heat  conduction 
Blood  perfusion 
Metabolism 
Circadian  rhythm 
Physical  activity 
Water  intake 
Convection 
Radiation 

Evaporative  cooling 
Temperature- 
dependent  changes  J 


Solution: 

Finite-element 

method 


Fig.  1.  Steps  in  the  development  of  the  computational  model.  The  preliminary 
anatomical  model  was  acquired  from  Duke  University.  The  anatomical  fea¬ 
tures  that  were  missing  from  the  preliminary  model  were  reconstructed  from 
stacks  of  magnetic  resonance  images  obtained  from  Brooks  Air  Force  Base. 
Preprocessing  of  the  geometry  for  subsequent  computations  using  finite- 
element  method  was  performed  in  Mimics  and  3-matic  (Materialise,  Plymouth, 
MI).  The  mathematical  model  was  formulated  by  incorporating  the  heat- 
transfer  mechanisms  shown  and  was  solved  in  COMSOL  Multiphysics 
(COMSOL,  Burlington,  MA).  MRI,  magnetic  resonance  imaging;  3-D,  three- 
dimensional. 


system  (6).  The  increases  in  skin  blood  flow  and  evaporation 
rate  due  to  perspiration/salivation  enhance  the  rate  of  heat  loss 
in  the  body,  as  it  attempts  to  cool  down.  Body  temperature  rises 
when  the  rate  of  heat  loss  is  unable  to  match  the  heat  gained 
due  to  thermal  stress  while  also  accounting  for  metabolic  heat 
generation  and  heat  transfer  away  from  the  organ.  To  accu¬ 
rately  quantify  the  temperature  distribution  in  the  body  both 
during  normothermic  and  thermal- stress  conditions,  key  ther¬ 
moregulatory  mechanisms  need  to  be  incorporated  into  the 
computational  model. 

Numerous  studies  have  focused  on  the  development  of 
whole  body  computational  models  of  thermoregulation  for 
animal  models  and  humans  (5,  27,  32,  35,  40,  44,  46,  71,  72). 
However,  most  of  these  studies  have  grossly  simplified  the 
anatomy  by,  for  example,  considering  a  one-dimensional  radial 
representation  or  assuming  regular  shapes,  such  as  cylinders 
and  spheres,  for  organs  and  body  contours  (27,  35,  40,  65). 
While  some  included  anatomically  realistic  description  of 
the  geometry,  the  emphasis  of  most  of  these  models  was  on  the 
short-term  evaluation  of  temperature  distribution  to  gauge  the 
effects  of  electromagnetic  irradiation  (5,  44,  71,  75,  76). 
Consequently,  owing  to  the  short  time  scales  (30-5,400  s)  of 
these  studies,  the  models  did  not  attempt  to  describe  longer 
periodic  phenomena  (>24  h),  such  as  circadian  variations,  or 
used  simplifying  assumptions,  such  as  constant  blood  temper¬ 
ature,  during  the  heating  process.  Furthermore,  some  models 
even  ignored  critical  heat  transfer  mechanisms,  such  as  con¬ 
vection  due  to  blood  perfusion  (71).  Heat  transfer  models  for 
other  applications,  such  as  burn  injury  and  thermal  therapy, 
have  focused  on  only  the  organ  or  region  of  interest  without 
including  the  whole  body  anatomy  (22,  23,  57,  90).  However, 
due  to  heat  transfer  by  cardiovascular  blood  flow  and  other 
mechanisms,  such  as  thermal  conduction,  thermal  responses  in 
the  various  organs  are  interdependent.  Moreover,  as  discussed 
above,  severe  heat  stress  may  lead  to  a  systemic  inflammatory 
response  and  varying  severity  of  dysfunctions  in  different 
organs.  Therefore,  to  investigate  the  progression  of  heat  stress 
and  elucidate  the  differential  response  in  different  organs,  it  is 
imperative  that  the  computational  model  characterizes  the 
whole  body  thermal  response  during  and  after  the  heat- stress 
period. 

In  this  work,  we  present  a  computational  approach  to  predict 
the  temporal  dynamics  and  spatial  temperature  distribution  for 
an  in  vivo  rat  model.  Our  computational  model  includes 
three-dimensional  (3-D)  geometry  obtained  from  medical  im¬ 
aging  to  accurately  characterize  the  rat  anatomy.  Our  mecha¬ 
nistic  computational  modeling  approach  incorporates  various 
mechanisms  of  heat  transfer  associated  with  thermal  stress, 
including  heat  conduction,  blood  perfusion,  metabolic  heat 
generation,  physical  activity,  water  intake,  evaporative  cooling, 
and  external  convection  and  radiation.  The  model  includes 
novel  formulations  for  computing  blood  temperature  as  well  as 
the  circadian  rhythm  in  body  temperature,  which  have  not  been 
considered  in  previous  efforts.  We  validated  the  predicted  core 
temperature  and  its  circadian  variation,  organ  volumes,  and 
organ  temperatures  using  in  vivo  animal  experiments  and 
literature  data.  The  model  provides  the  capability  to  predict  the 
temperature  distribution  within  each  vital  organ  and  potentially 
establish  correlations  with  impending  inflammatory  responses 
and  differential  organ  dysfunction.  This,  in  turn,  will  facilitate 
the  identification  of  organ-specific  risks  due  to  exposure  to 
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thermal  stress,  such  as  acute  liver  failure  and  brain  inflamma¬ 
tion  (13,  31).  In  addition,  it  would  potentially  aid  in  the 
development  of  improved  clinical  strategies  for  prevention  and 
management  of  heat  stress,  such  as  design  of  localized  cooling 
strategies  (47). 

Glossary 

A  Surface  area  (m2) 

AL  Activity  level  (dimensionless) 

ac  Relative  activity  level  (dimensionless) 

acmajc  Maximum  relative  activity  level  (dimensionless) 
amet  Amplitude  of  circadian  variation  in  metabolic 
heat  generation  (dimensionless) 
aper  Amplitude  of  circadian  variation  in  blood  per¬ 
fusion  (dimensionless) 
cp  Specific  heat  capacity  (J-kg_1-°C_1) 

cPyb  Specific  heat  capacity  of  blood  (J-kg_1-°C_1) 

cp, water  Specific  heat  capacity  of  water  (J-kg_1-°C_1) 

Cmet  Circadian  function  for  metabolic  heat  genera¬ 
tion  (dimensionless) 

Cper  Circadian  function  for  blood  perfusion  (dimen¬ 
sionless) 

d  Characteristic  diameter  (m) 

hc  Convective  heat  transfer  coefficient  (W-m_2oC_1) 

k  Thermal  conductivity  (W-m_1-°c_1) 

Lwater  Latent  heat  of  vaporization  of  water  (J/kg) 

M  Body  mass  (kg) 

Mwater  Experimentally  measured  water  intake  (kg) 
riisaiiva  Rate  of  saliva  loss  (kg/s) 

n  Direction  normal  to  the  skin  surface  (dimen¬ 
sionless) 

Qm  Metabolic  heat  generation  rate  (W/m3) 

Qmo  Basal  metabolic  heat  generation  rate  (W/m3) 
Qioss  Volumetric  heat  loss  due  to  water  intake  (W/m3) 
T  Body  temperature  at  a  given  x,  y,  z  location  (°C) 
t  Time  (s) 

To  Reference  temperature  at  steady  state  (°C) 

Ta  Arterial  blood  temperature  (°C) 

Tamb  Ambient  temperature  (°C) 

Tc  Core  temperature  (°C) 

Tl  Body-temperature  limit  corresponding  to  base¬ 

line  blood  perfusion  (°C) 

Tskin  Skin  surface  temperature  (°C) 

Tu  Body-temperature  limit  corresponding  to  ele¬ 
vated  or  depressed  blood  perfusion  (°C) 

T water  Water  temperature  (°C) 

u  Air  flow  velocity  (m/s) 

V  Organ  volume  (m3) 

V£q  Basal  blood  perfusion  rate  (m3-s_1-m-3) 

Vf  Blood  perfusion  rate  (m3-s_1-m_3) 

VlHS  Blood  perfusion  rate  under  heat-stress  condi¬ 

tions  (m3-s_1-m_3) 
y stomach  Volume  of  stomach  (m3) 

a  Coefficient  for  linear  increase  in  perfusion  with 
temperature  (dimensionless) 

/3  Coefficient  for  linear  decrease  in  perfusion  with 
temperature  (dimensionless) 

A t  Time  step  size  (s) 

A twater  Water  intake  time  interval  (s) 

s  Emissivity  (dimensionless) 


7)  Mechanical  efficiency  in  metabolic  heat  gener¬ 
ation 

/x  Dynamic  viscosity  (Pa-s) 

p  Density  (kg/m3) 

pb  Density  of  blood  (kg/m3) 

cr  Stefan-Boltzmann  constant  (5.67  X  10-8 
W-m_2oC-4) 

t  Period  of  the  circadian  oscillator  (s) 
cpmet  Circadian  phase  for  metabolic  heat  generation 
(dimensionless) 

(Pper  Circadian  phase  for  blood  perfusion  (dimensionless) 

MATERIALS  AND  METHODS 

The  computational  model  was  formulated  by  including  a  3-D 
representation  of  the  rat  geometry  obtained  from  magnetic  resonance 
imaging  (MRI)  and  describing  the  dominant  modes  of  heat  transfer 
during  thermoregulation.  Simulations  were  performed  for  control  and 
heat- stress  conditions  that  mimicked  the  experimental  conditions.  The 
details  of  model  development  steps  are  illustrated  in  Fig.  1  and 
discussed  below  along  with  the  experimental  protocol. 

Experimental  Protocol 

Simulations  were  set  up  based  on  in  vivo  rat  experiments  per¬ 
formed  at  the  U.S.  Army  Research  Institute  of  Environmental  Medi¬ 
cine  (USARIEM).  The  Institutional  Animal  Care  and  Use  Committee 
approved  all  experimental  procedures,  which  were  performed  in 
accordance  with  the  American  Physiological  Society’s  Guiding  Prin¬ 
ciples  for  Research  Involving  Animals  and  adhered  to  the  National 
Institutes  of  Health  Guide  for  Care  and  Use  of  Laboratory  Animals. 

Animals.  A  power  analysis,  performed  using  core  temperature  data 
from  a  previous  rat  heat- stress  study,  indicated  that  n  =  15  control  rats 
and  n  =  15  heat-stressed  rats  were  needed  to  achieve  statistical 
significance  for  observing  mean  differences  of  0.5 °C  between  the  two 
groups  during  recovery  after  heat  stress.  However,  in  this  study,  we 
leveraged  data  for  model  development  and  validation  from  a  larger 
study  ( n  =  80)  designed  to  establish  a  plasma  and  tissue  bank  for 
control  (n  =  30)  and  surviving  (>1  day)  heat-stressed  (n  =  50) 
animals.  Male  Fischer  344  (F344;  n  =  80;  Charles  River  Laboratories, 
Stone  Ridge,  NY)  rats  weighing  234-336  g  (—2-3  mo  old)  were  used. 
Animals  were  obtained  in  five  separate  shipments,  with  rats  randomly 
assigned  to  heat- stressed  and  control  groups.  Rats  were  individually 
housed  in  Nalgene  polycarbonate  cages  (10.5  X  19  X  8  in.,  Ancare, 
Bellmore,  NY)  with  HEPA-filter  cage  tops  and  ALPHA-dri/Cob  blend 
bedding  (PharmServ,  Framingham,  MA).  Environmental  enrichment 
consisted  of  a  rat  Igloo  (Nalgene  Nunc,  Rochester,  NY),  a  small 
plastic  ball,  a  Nylabone  for  chewing,  and  a  Maplewood  product 
containing  a  food  treat  to  encourage  foraging  (W0002,  Bio-Serv, 
Frenchtown,  NJ).  Rats  were  housed  under  standard  laboratory  condi¬ 
tions  (22.0  ±  0.2°C,  12:12-h  light-dark  cycle,  lights  on  at  0600)  in 
an  Association  for  Assessment  and  Accreditation  of  Laboratory  An¬ 
imal  Care-accredited  facility.  A  relatively  cool  housing  temperature 
(22.0  ±  0.2°C)  was  chosen  for  these  experiments  to  ensure  rapid 
cooling  and  survival  during  heat-stress  recovery.  Rodent  laboratory 
chow  (LM-485,  Harlan  Teklad,  Madison,  WI)  and  water  were  pro¬ 
vided  ad  libitum  except  during  experimental  procedures.  Water  in¬ 
gested  by  the  animals  was  measured. 

Radiotelemetry  measurements.  Rats  were  implanted  with  TL1 1M2- 
C50-PXT  PhysioTel  Multiplus  Transmitters  (Data  Sciences,  St.  Paul, 
MN)  10  days  before  arrival  at  USARIEM.  The  temperature  sensors 
were  implanted  in  the  intraperitoneal  cavity  and  were  sutured  to  the 
peritoneal  wall.  The  transmitters  measured  core  temperature  (77; 
±0.25°C)  and  did  not  produce  thermal  energy  on  implantation.  Each 
transmitter  emitted  a  unique  frequency  that  was  proportional  to  Tc  and 
was  received  by  an  antenna  under  the  cage  of  each  animal.  All 
radiotelemetry  variables  were  monitored  at  1-min  intervals  in  con- 
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scious,  free-ranging  animals  that  maintained  normal  heat  dissipation 
mechanisms.  Implanted  transmitters  weighed  ~11  g  or  —3-4%  of  the 
experimental  body  weight  of  the  animal.  Animals  were  housed  for  >  1 
wk  before  experimental  protocols  to  acclimatize  to  the  environmental 
conditions  and  ensure  a  robust  circadian  Tc  rhythm  before  experimen¬ 
tation.  We  calculated  the  average  daytime  Tc  values  (0600-1759)  for 
all  rats  for  2  days  prior  to  experimentation  (after  the  acclimatization 
period  to  our  facility)  and  used  this  value  (37.3 °C)  as  the  criterion  for 
baseline  temperature  at  which  heat  stress  could  be  initiated  for  each 
animal. 

Heating  protocol.  All  experiments  were  conducted  in  conscious, 
free-moving  animals.  Randomly  assigned  rats,  in  their  original  cage 
with  filter  tops  removed,  were  placed  in  a  floor  standing  incubator 
(Thermo  Scientific,  Ashville,  NC)  set  at  room  temperature  (22.0  ± 
0.2°C)  24  h  before  the  initiation  of  heat-stress  experiments.  Non- 
heated  rats  were  not  introduced  to  the  incubator  environment  and 
remained  in  their  original  cage  and  rack  location  throughout  the 
heating  protocols.  Heat-stress  experimentation  was  initiated  between 
0800  and  1000,  after  the  Tc  of  each  rat  reached  values  below  37.3°C 
(the  12-h  daytime  average  Tc  of  all  animals  as  discussed  above).  The 
heating  experiments  were  initiated  by  weighing  rats  in  the  heat- 
stressed  ( n  =  50)  and  control  (; n  =  30)  groups  followed  by  removal  of 
food  and  water.  Rats  in  the  heat- stressed  group  were  placed  in  the 
incubator  in  their  original  cage,  and  the  incubator  temperature  was 
increased  to  37.0  ±  0.2°C.  Control  animals  were  placed  in  their 
original  location  at  the  normal  housing  temperature  of  22.0  ±  0.2°C. 
Experimental  animals  were  heated  until  a  Tc  of  41.8°C  ( TcMax )  was 
reached.  After  attainment  of  TcMax,  rats  were  removed  from  the 
incubator,  weighed,  placed  in  a  new  cage,  and  returned  to  normal 
housing  temperature  (Ta  =  22.0  ±  0.2°C).  Time-matched  control  rats 
were  weighed,  provided  a  new  cage,  and  placed  back  at  their  original 
cage  location.  Control  and  experimental  animals  were  provided  food 
and  water  ad  libitum  throughout  recovery.  Body  weights  were  ob¬ 
tained  on  a  top-loading  balance  with  an  accuracy  of  ±0.1  g. 

Rat  Anatomical  Model 

The  3-D  anatomical  model  of  a  rat  was  acquired  from  Duke 
University  (49).  The  model  was  developed  by  segmentation  and  3-D 
reconstruction  of  medical  images  obtained  from  MRI  based  on  non- 
uniform  rational  B- splines  and  includes  a  realistic  description  of  the 
organs  and  skeletal  system  of  the  rat  (Fig.  1).  Missing  anatomical 
features,  such  as  the  tail,  were  added  to  the  model  by  3-D  reconstruc¬ 
tion  of  MR  images  (81  images,  3-mm  thick)  obtained  from  Brooks  Air 
Force  Base  (58)  using  the  Mimics  (Materialise,  Plymouth,  MI)  med¬ 
ical  image  processing  software  package.  The  native  format  of  the 
anatomical  model  was  the  standard  stereolithography  file  format, 
which  describes  the  surface  of  geometry  using  a  mesh  of  triangles. 
The  primary  objective  of  the  original  anatomical  model  acquired  from 
Duke  University  was  visualization.  Therefore,  while  the  model  accu¬ 
rately  described  the  rat  anatomy  in  3-D,  it  had  artifacts  that  made  it 
incompatible  with  the  development  of  a  computational  model  using 
finite-element  method  (FEM),  especially  during  the  volume  mesh 
generation  process.  Consequently,  the  anatomical  model  was  im¬ 
ported  into  a  FEM  preprocessor  and  computer-aided  design  (CAD) 
software  package,  3-matic  (Materialise),  and  artifacts,  such  as  inter¬ 
secting  organs  and  bones,  were  removed.  The  anatomical  model  had 
other  artifacts,  such  as  presence  of  a  large  number  of  bad-quality  and 
inverted  triangles,  and  local  distortions  in  the  geometry  that  made  it 
incompatible  with  the  development  of  a  volume  mesh  for  the  com¬ 
putational  model.  These  artifacts  were  removed  by  local  smoothing 
and  meshing  of  the  geometry  in  3-matic.  The  improved  surface  mesh 
was  then  used  to  create  the  volume  mesh.  The  mesh  consisted  of 
451,762  tetrahedral  volume  elements,  as  shown  in  Fig.  2.  The  pro¬ 
cessed  whole  body  anatomical  model  included  essential  organs,  such 
as  the  heart,  liver,  kidneys,  brain,  lungs,  stomach,  and  intestines,  as 
well  as  the  skeletal  system.  The  anatomical  model  with  the  mesh  was 
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Lungs  Liver  Intestines 


Fig.  2.  Geometry  (A)  and  mesh  (. B )  of  the  3-D  rat  anatomical  model.  The 
geometry  includes  60  volumes  and  70  boundaries  representing  the  organs, 
interfaces,  and  external  surfaces.  The  various  organs  included  in  the  model  are 
shown  in  different  colors.  The  mesh  consists  of  451,762  elements. 

imported  into  a  FEM-based  software  package,  COMSOL  Multiphys¬ 
ics  (COMSOL,  Burlington,  MA),  to  develop  the  computational 
model. 

Computational  Model 

The  physiological  model  for  thermoregulation  in  rats  during  heat 
stress  was  formulated  based  on  the  Pennes  bioheat  transfer  equation 
(64,  85).  The  Pennes  equation  was  chosen  because  it  has  been 
successfully  employed  in  numerous  whole  body  human  and  animal 
heat  transfer  models  (27,  40,  44).  The  equation  describes  the  transient 
change  in  temperature  at  any  location  inside  the  rat  and  accounts  for 
heat  conduction,  heat  convection  due  to  blood  perfusion  in  the  tissues, 
and  metabolic  heat  generation.  The  Pennes  equation  approximates  the 
heat  transfer  due  to  blood  perfusion  in  the  body  organs  as  a  source  (or 
sink)  term  instead  of  explicitly  modeling  blood  flow  and  heat  ex¬ 
change  in  the  blood  vessels,  as  follows: 

pUU  V-(feVr)  +  Pbcp,bVb(Ta-T)  +  Qm  ( 1 ) 

o  t  ^  ^  ^  J  ^ 

heat  conduction  heat  convection  metabolic 

due  to  blood  perfusion  heat  generation 

where  p,  cp,  and  k  denote  tissue  density,  specific  heat  capacity,  and 
thermal  conductivity,  respectively;  pb  and  cp>h  represent  the  density  and 
specific  heat  capacity  of  the  blood,  respectively;  Vh  denotes  the  blood 
perfusion  rate  in  the  tissue,  and  Ta  represents  the  arterial  blood  temperature. 

Circadian  variations.  Because  one  of  the  aims  of  the  model  was  to 
predict  the  temperature  distribution  in  heat-stressed  rats  during  the 
course  of  a  day,  it  was  critical  to  incorporate  the  effects  of  circadian 
rhythmicity.  Both  blood  perfusion  and  metabolic  heat  generation  rates 
demonstrate  circadian  variations  in  rats,  as  reported  in  the  literature 
(21,  68,  73,  82).  The  circadian  process  for  blood  perfusion  was 
modeled  by  implementing  the  following  sinusoidal  function: 

/  \  .  /  2'nT  \ 

Cperyt)  ^ per  SUll  ^  T  (fper  J  (2) 

where  aper  denotes  the  amplitude  of  circadian  variation  in  blood 
perfusion,  c pper  represents  the  circadian  phase,  and  t  denotes  the 
period  of  the  circadian  oscillator  (24  h).  The  amplitude  of  circadian 
variation  in  blood  perfusion  and  the  circadian  phase  for  each  organ 
were  determined  by  fitting  experimental  data  on  blood  perfusion  (21) 
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with  Eq.  2.  The  circadian  rhythmicity  in  metabolic  heat  generation 
was  modeled  similarly,  as  follows: 


Cmetif) 


/27Tt  \ 

^ met  SHI  I  ^  +  (Pmet  J 


(3) 


where  amet  denotes  the  amplitude  of  circadian  variation  in  metabolic 
heat  generation,  c pmet  represents  the  circadian  phase,  and  t  denotes  the 
period  of  the  circadian  oscillator  (24  h).  Experimental  data  on  meta¬ 
bolic  heat  generation  (68)  were  used  to  obtain  the  circadian  variation 
and  the  circadian  phase  in  Eq.  3. 

Blood  flow  model.  The  Pennes  bioheat  equation  (Eq.  1 )  includes  the 
blood  perfusion  term  as  a  source  (or  sink)  term  in  the  standard  heat 
equation,  with  the  rate  of  heat  transfer  between  blood  and  tissue 
depending  on  the  blood  perfusion  rate  and  the  difference  between 
arterial  blood  and  tissue  temperatures.  In  this  formulation,  it  is 
assumed  that  the  heat  exchange  between  blood  and  tissue  takes  place 
in  the  capillary  beds  and  that  the  blood  equilibrates  with  the  tissue  so 
that  the  temperature  of  the  blood  flowing  out  of  the  veins  is  equal  to 
the  local  tissue  temperature  (60).  Pennes  (64)  did  not  incorporate  the 
subsequent  mixing  of  the  venous  blood  from  various  organs  and 
tissues  in  his  original  work  because  his  model  included  the  forearm 
only.  However,  for  a  whole  body  model,  the  arterial  blood  tempera¬ 
ture  would  depend  on  the  mixing  of  the  venous  blood  from  different 
tissues  before  it  is  pumped  back  by  the  heart  to  the  pulmonary  and 
systemic  circulations.  Incorporating  the  arterial  blood  temperature 
(Ta)  as  a  constant  value,  as  implemented  in  the  original  formulation  by 
Pennes  (64)  and  subsequent  work  (23,  57,  90),  can  only  be  justified 
when  modeling  a  local  subsystem  of  the  body  (e.g.,  an  arm)  or,  if 
modeling  the  entire  body,  for  only  a  short  duration  of  heating  or 
cooling,  when  the  global  arterial  temperature  is  not  expected  to 
change  significantly.  For  whole  body  thermoregulation  models  that 
aim  to  investigate  long-heating  periods,  as  in  the  proposed  work,  the 
assumption  of  a  constant  blood  temperature  is  not  valid.  Therefore,  we 
used  a  novel  mixing  formulation  for  estimating  the  arterial  tempera¬ 
ture  at  any  time  t ,  as  follows: 


2  /  VL-7TA'  dv, 


where  the  subscript  i  refers  to  the  various  organs  included  in  the 
model  and  the  time  interval  (At)  is  the  time  step  size  for  the 
computational  model.  Equation  4  calculates  the  volumetric  flow- 
weighted  average  temperature  for  the  entire  body  based  on  the 
temperature  distribution  in  the  body  organs  at  the  last  computed  time 
step. 

Temperature  dependence  of  blood  flow.  Blood  perfusion  has  been 
reported  to  be  dependent  on  temperature,  with  the  rate  doubling  for 
each  10°C  rise  in  temperature  from  the  baseline  (35).  This  relationship 
was  implemented  in  the  model  using  the  following  equation: 


VI  =  VL  2V 


(5) 


Equation  5  was  successfully  used  to  model  human  thermoregulation 
(35),  and  a  modified  version  of  the  equation,  also  derived  from  human 
models  (5),  was  used  for  modeling  rats  (44). 

Additionally,  Kregel  at  al.  (51)  studied  the  circulatory  dysfunction 
and  observed  a  differential  thermoregulatory  response  in  splanchnic 
and  peripheral  blood  flows  during  heat  stress.  The  splanchnic  blood 
flow  was  observed  to  decrease  during  heat  stress  as  opposed  to  an 
increase  in  the  peripheral  blood  flow.  In  the  absence  of  detailed 
measurements  of  changes  in  blood  flow  in  the  different  organs  due  to 
heat  stress,  we  implemented  a  linear  dependence  of  blood  flows  on 
organ  temperatures  based  on  the  works  of  Kregel  et  al.  (51)  and  Hirata 
et  al.  (44),  as  follows: 
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In  this  formulation,  blood  perfusion  remained  at  a  baseline  level 
below  a  certain  temperature  (Tf)\  linearly  increased  (in  the  outer  skin 
layer  and  tail)  or  decreased  (in  the  internal  organs)  with  temperature 
for  a  fixed  temperature  range;  and  remained  at  an  elevated  or  de¬ 
pressed  level  after  the  upper  bound  temperature  (7V)  was  reached. 
Equation  6  has  been  previously  used  for  modeling  both  rat  and  human 
thermoregulation  (5,  44,  45). 

Metabolic  heat  generation  rate.  Body  organs  are  known  to  have 
significantly  different  metabolic  activity  levels.  For  example,  or¬ 
gans  such  as  the  liver,  brain,  kidneys,  and  heart  have  been  observed 
to  be  highly  metabolically  active  compared  with  other  organs  and 
tissues  (25,  36,  81).  To  model  the  differential  thermal  response  in 
organs,  it  is  imperative  to  accurately  account  for  the  metabolic  heat 
generation  term  of  each  organ.  Wang  et  al.  (81)  used  Kleiber’s  law 
(50),  which  relates  the  total  basal  metabolic  heat  generation  rate 
with  body  mass,  to  derive  relationships  between  the  basal  meta¬ 
bolic  rates  of  different  organs  and  body  mass.  Metabolic  heat 
generation  rates  in  our  model  were  calculated  using  these  relation¬ 
ships  and  are  listed  in  Table  1. 

Temperature  dependence  of  the  metabolic  heat  generation  rate. 
The  metabolic  heat  generation  rate  in  any  organ  is  a  function  of 
biochemical  reactions  occurring  in  the  cells.  This  reaction  rate  doubles 
for  each  10°C  rise  in  temperature  from  the  baseline.  Accordingly,  the 
temperature  dependence  of  the  metabolic  heat  generation  rate  was 
given  by  the  following  (84): 

(hl°) 

QZsan(T)  =  Qzr^  w  1  (7) 

Heat  generation  due  to  physical  activity.  Physical  activity  of  the 
rats  inside  an  incubator  was  measured  by  recording  the  activity 
count  [ac(0]  using  radiotelemetry  (Fig.  3).  Activity  count  is  a 
relative  measure  of  the  distance  and  speed  of  movement  of  the 
animal.  Heat  generated  due  to  physical  activity  was  modeled  by 
incorporating  the  measured  activity  data  in  the  mathematical  for¬ 
mulation.  However,  since  activity  count  provides  only  a  relative 
measure  of  physical  activity,  the  following  assumption  was  em¬ 
ployed  to  quantify  the  effect  of  physical  activity:  because  rats 
inside  an  incubator  were  restricted  to  a  light  activity  level,  the 
maximum  measured  activity  count  was  assumed  to  be  correlated  to 
the  maximum  metabolic  heat  generation  rate  associated  with  light 
activity,  as  obtained  from  the  literature  (1,  8,  9).  Subsequently, 
accounting  for  the  rate  of  thermal  energy  generation  converted  to 
mechanical  power,  the  following  relationship  was  used  to  estimate 
the  heat  generation  rate  due  to  activity  in  the  model  (27): 

,  .  ac (t) 

Qm,activity  =  ALgm(l  -  T }) -  (8) 

^max 


The  activity  level  (AL)  for  light  activity  and  the  corresponding 
mechanical  efficiency  in  the  conversion  of  metabolic  heat  into  me¬ 
chanical  energy  (p)  were  obtained  from  the  literature  and  are  listed  in 
Table  2.  The  average  measured  activity  counts  [ac(0]  for  the  control 
(n  =  30)  and  heat-stressed  (n  =  50)  rats  were  calculated  (Fig.  3)  and 
input  to  the  model  as  a  function  of  time. 
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Tissue 

Properties,  Units 

Reference(s) 

p,  kg/m3 

k,  W-m_1-°C_1 

cp,  J-kg  1-°C  1 

Vbov,  m3-s  x-m  3 

aPer 

Qmo,  W/m3 

Stomach 

1,040 

0.518 

3,500 

0.007 

0.420 

1,837 

21,  71,  81 

Kidneys 

1,040 

0.502 

3,638 

0.107 

0.240 

38,128 

21,  71,  81 

Muscle 

1,047 

0.518 

3,500 

0.001 

0.170 

1,849 

19,  75,  81 

Bone 

1,600 

0.300 

2,000 

0.030 

0.190 

2,826 

21,  75,  81 

Brain 

1,036 

0.528 

3,854 

0.006 

0.210 

26,351 

21,  75,  81 

Fat 

970 

0.204 

4,483 

0.001 

0.570 

1,713 

19,  75,  81 

Heart 

1,030 

0.526 

3,329 

0.049 

0.200 

51,041 

21,  75,  81 

Liver 

1,030 

0.511 

3,871 

0.022 

0.200 

46,658 

19,  75,  81 

Lungs 

1,040 

0.252 

1,935 

0.253 

0.200 

1,837 

63,  71,  81 

Tendon 

1,120 

0.500 

3,118 

0.002 

0.200 

1,978 

26,  43,  52,  81 

Skin 

1,100 

0.370 

3,700 

0.002 

0.280 

1,943 

19,  75,  81 

Intestines 

1,042 

0.518 

3,500 

0.023 

0.200 

1,841 

21,  75,  81 

Blood 

1,057 

3,600 

18 

See  Glossary  for  definitions. 


Heat  exchange  due  to  water  intake.  Rats  were  freely  allowed  to  drink 
water  kept  at  room  temperature  at  all  times  except  during  heat  stress. 
Accordingly,  the  heat  loss  due  to  this  water  intake  was  modeled  as  a  sink 
term  in  the  Pennes  bioheat  equation  for  the  stomach,  as  follows: 
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Fig.  3.  The  relative  motor  activities  of  control  rats  (n  =  30)  placed  at  an 
ambient  temperature  of  22°C  (A)  and  heat-stressed  rats  (n  =  50)  in  an 
incubator  maintained  at  37°C  (. B )  for  163  min.  The  solid  lines  represent  the 
average  of  the  measured  activity  values  and  the  gray  bars  indicate  one  standard 
deviation.  Motor  activity  data  were  used  as  input  parameters  to  the  model  to 
modulate  the  metabolic  heat  generation  rate.  Light  and  dark  cycles  are 
indicated  by  white  and  black  bars,  respectively. 
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where  mwater  was  the  water  ingested  by  the  animal  as  measured  in  the 
experiment. 

External  convection.  Under  control  conditions,  rats  were  kept  at 
normal  room  temperature  and,  therefore,  their  primary  mechanism  of 
heat  transfer  with  the  surroundings  was  convection.  Similarly,  rats 
were  heated  by  convection  due  to  laminar  airflow  inside  an  incubator 
set  at  37°C  under  heat-stress  conditions.  Heat  transfer  due  to  convec¬ 
tion  was  included  in  the  model  as  a  boundary  condition  applied  to  the 
surface  of  the  rat’s  skin  as  follows: 


d  T  ki 

-k— —  =  hc(T-Tamb) 
on 


(10) 


To  emulate  the  experimental  conditions  for  control  and  heat- stressed 
animals,  different  ambient  temperatures  ( Tamb )  and  heat  transfer 
coefficients  (hc)  were  used,  as  listed  in  Table  2.  Heat  stress  in  the 
computational  model  was  prescribed  by  applying  the  convective 
boundary  condition  for  the  same  application  time  as  per  the  experi¬ 
mental  protocol.  Heat  transfer  coefficients  for  the  control  and  heat- 
stress  cases  were  obtained  using  previously  developed  correlations 
available  in  standard  heat  transfer  literature  (48).  The  corresponding 
heat  transfer  coefficients  for  the  head  and  trunk,  forelimbs,  hindlimbs, 
and  tail  of  the  rat  were  obtained  using  the  following  equation  (48): 

_  0.683  p0A66u0A66c0333k0.661 
h c  ~  jO.534  ^0.133 

The  characteristic  diameters  (d)  of  the  various  body  parts  were 
calculated  based  on  their  volumes,  and  the  air  properties  (p,  cp,  k,  /x) 
were  obtained  from  the  literature  (18)  for  the  two  ambient  tempera¬ 
tures  corresponding  to  control  and  heat-stress  conditions. 

Radiation.  In  addition  to  convection,  the  rat’s  body  surface  ex¬ 
changed  heat  with  the  surroundings  due  to  radiation.  Radiative  heat 
transfer  was  modeled  by  applying  the  following  boundary  condition 
on  the  rat’s  body  surface: 

-kd-^  =  s*(T*-Tlb)  (12) 

on 

Evaporative  cooling.  Due  to  the  absence  of  sweat  glands,  rats  do 
not  perspire  (28).  The  ability  to  spread  saliva  on  the  surface  of  their 
body  and  the  resulting  heat  transfer  due  to  evaporation  is,  therefore,  a 
critical  heat  regulation  mechanism  for  rats.  Damas  (17)  reported  that 
evaporative  water  loss  corresponds  to  <10%  of  the  total  body  weight 
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Table  2.  Input  and  model  parameters  in  addition  to  the 
thermal  and  blood  flow  model  properties  listed  in  Table  1 


Parameter,  Unit 

Value 

Reference(s) 

Input  parameters 

Ambient  temperature  (control  condition) 

(Tomb),  °C 

22 

Experiment 

Ambient  temperature  (during  heat  stress) 

(Tamb),  °C 

37 

Experiment 

Water  intake  (mwater),  kg 

Experiment 

Control  conditions 

16.47  X  10“3 

During  heat  stress 

29.38  X  10“3 

Rate  of  saliva  loss  (msaiiva),  kg/s 

Experiment 

Control  conditions 

4.41  X  10-8 

During  heat  stress 

2.34  X  1(T7 

Heating  time,  s 

9780 

Experiment 

Water  temperature  (Twater),  °C 

22 

Experiment 

Model  parameters 

Reference  temperature  (T0),  °C 

39 

44 

Specific  heat  capacity  of  water  ( cp,water ), 

4,180 

18 

Convective  heat  transfer  coefficient  (hc). 

W-m_2-°C_1 

38 

Control  conditions 

Head  and  trunk 

4.0 

Forelimbs 

6.2 

Hindlimbs 

4.6 

Tail 

9.0 

During  heat  stress 

Head  and  trunk 

16.2 

Forelimbs 

24.9 

Hindlimbs 

18.6 

Tail 

35.8 

Emissivity  (e),  dimensionless 

0.8 

76 

Latent  heat  of  vaporization  of  water 

( Lwater )?  J/kg 

2.26  X  106 

18 

Activity  level  for  light  activity  (AL), 

dimensionless 

1.6 

1 

Mechanical  efficiency  at  light  activity 

(rj),  dimensionless 

0 

27 

Amplitude  of  circadian  variation  in 

metabolism  (amet) 

0.2 

68 

Period  of  the  circadian  oscillator  (t),  s 

86,400 

21,  68 

Circadian  phase  (<p),  dimensionless 

0 

21,  68 

These  input  and  model  parameters  were  either  obtained  from  the  literature 
or  were  based  on  the  experimental  conditions,  as  noted. 


loss  in  rats  at  an  ambient  temperature  of  22°C  and  —70%  for  rats  heat 
stressed  at  37°C,  with  urination,  defecation,  and  oxidative  metabolism 
accounting  for  the  remaining  weight  loss.  Accordingly,  total  body 
fluid  loss  due  to  salivation  was  estimated  based  on  the  total  body 
weight  loss  of  the  rats  measured  during  the  experiment.  The  rate  of 
saliva  loss  was  assumed  to  vary  linearly  with  temperature,  as  reported 
in  the  literature  (29,  74).  Finally,  heat  loss  per  unit  area  due  to 
salivation  was  included  in  the  model  as  the  following  boundary 
condition: 

,  ^  Tsjcin  ifl salivaflwater 

—  k - =  -  (13) 

dn  Arat 

Solution  details.  The  computational  model,  including  the  3-D 
anatomical  model,  governing  Eq.  7,  auxiliary  equations  (Eqs.  2-9), 
and  boundary  conditions  (Eqs.  10-13),  were  set  up  and  solved  using 
the  FEM  solver,  COMSOL  Multiphysics.  Input  and  model  parameters 
for  the  simulations  are  listed  in  Tables  1  and  2.  Steady-state  simula¬ 
tions  were  performed  to  initialize  the  temperature  distribution  for 
control  and  heat- stressed  rats  at  the  beginning  of  the  heating  protocol. 
For  the  steady-state  simulation,  boundary  conditions  corresponding  to 
normal  room  temperature  were  applied.  After  the  temperature  distri¬ 


bution  in  the  computational  domain  (i.e.,  the  rat  geometry)  was 
initialized,  time-dependent  simulations  were  carried  out  to  predict 
temperature  distributions  over  a  25 -h  time  period.  Heating  conditions 
corresponding  to  an  ambient  temperature  of  37°C  for  163  min  were 
applied  in  the  model  for  simulating  heat-stressed  rats  based  on  the 
experimental  protocol. 

RESULTS 

Model  Validation 

The  first  step  in  model  validation  was  to  compare  the  model- 
computed  variables  with  the  corresponding  data  from  the  litera¬ 
ture,  as  shown  in  Table  3.  However,  due  to  the  unavailability  of 
complete  anatomical  and  physiological  data  for  F344  rats,  data 
from  various  rat  strains  were  used  for  comparisons.  Organ  vol¬ 
umes,  which  were  computed  by  volume  integration  in  the  FEM 
solver,  were  found  to  be  consistent  with  those  reported  in  the 
literature  (19,  25,  30,  62,  63,  86),  establishing  the  anatomical 
accuracy  of  our  model.  Subsequently,  physiological  variables, 
such  as  core  and  organ  temperatures  and  total  metabolism  and 
circadian  variations,  predicted  by  our  model  were  compared  with 
those  reported  in  the  literature,  and  a  fairly  good  match  was 
observed  (Table  3).  The  most  striking  features  of  the  comparisons 
were  the  accurate  predictions  of  relative  temperatures  in  the 


Table  3.  Comparison  of  our  computational  model  with 
anatomical  and  physiological  data  for  normal  control  rats 
reported  in  the  literature 


Variable 

Model 

Computed 

Experimental 

Data 

Reference(s) 

Surface  area/body  mass, 

m2/kg 

0.12 

0.13 

33 

Organ  volumes,  X 10  6  m3 

Stomach 

4.70 

1.06-4.58 

Kidneys 

2.62 

2.00-3.70 

Brain 

2.05 

1.20-2.07 

19,  25,  30,  62,  63,  86 

Heart 

2.23 

1.20-2.35 

Liver 

1.04 

1.07-1.96 

Lungs 

2.93 

1.54-2.10 

Intestines 

11.50 

10.60-15.80 

Core  temperature  (Tc),  °C 

37.20 

37.0-37.4 

33,  68 

Circadian  variation  in  Tc,  °C 

0.80 

0.34-1.20 

33,  59,  68 

Metabolism,  W 

1.93 

1.30-1.99 

33,  68 

Metabolism,  W/kg 

6.20 

4.5-7 .4 

33 

Circadian  variation  in 

metabolism,  W 

0.77 

0.70 

68 

Organ  temperatures 

compared  with  blood 
temperature,  °C 

Liver 

0.40 

0.2-0. 3 

14,  15 

Kidneys 

0.20 

0.2-0. 3 

14,  15 

Intestines 

0.10 

0.1-0.13 

14,  15 

Hindlimbs 

-0.50 

-0.1  to 
-0.3 

14,  15 

Skin 

-0.70 

-0.6  to 
-1.2 

14,  15 

Tail 

-12.80 

-10.2  to 
-12.0 

34,  70 

Surface  area  and  organ  volumes  were  determined  by  computing  the  surface 
and  volume  integrals  of  the  geometry  in  the  finite-element  method  solver.  Core 
temperature  was  calculated  using  Eq.  1.  Total  metabolism  was  determined  by 
integrating  the  metabolic  heat  generation  rate  at  each  point  in  the  geometry, 
given  by  the  basal  metabolic  heat  generation  rate  (listed  in  Table  1)  and  Eqs. 
3,  7,  and  8.  Organ  temperatures  were  the  mean  temperatures  obtained  by 
calculating  the  volume-weighted  average  of  the  temperature  field  in  each  organ 
given  by  Eq.  1. 
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organs,  which  were  consistent  with  published  results  (14,  15) 
under  normothermic  conditions:  the  liver,  kidneys,  and  intestines 
were  predicted  to  be  at  a  higher  temperature  compared  with  the 
blood,  whereas  the  hindlimbs,  skin,  and  tail  were  cooler. 

For  a  quantitative  validation  of  the  model  dynamics,  pre¬ 
dicted  core  temperatures  as  a  function  of  time  of  day  were 
compared  with  those  obtained  experimentally  using  radiote¬ 
lemetry  measurements,  as  shown  in  Fig.  4.  For  the  control  rats, 
Fig.  4 A  shows  that  the  model  predictions  were  within  1  SD  of 
the  measured  data  ( n  =  30)  and  accurately  captured  the 
amplitude  of  the  circadian  rhythmicity  in  core  temperatures 
during  the  course  of  a  day-night  cycle.  Extreme  variations  in 
core  temperatures  (— 1.0°C)  in  control  rats  close  to  the  start  of 
the  experiments  (approximately  between  0800  and  1100), 
which  could  not  be  captured  by  the  model,  were  attributed  to 
the  external  disturbance  to  the  animals  due  to  weighing  and 
removal  of  food  and  water  at  the  start  of  the  heating  regime  and 
later  during  refeeding  after  the  heating  period  (to  mimic  the 
disturbances  experienced  by  heat- stressed  animals).  The  pre¬ 
dicted  core  temperature  decreased  as  the  day  progressed  and 
then  peaked  during  the  night  with  a  peak-to-peak  amplitude  of 
0.8°C.  However,  the  model  was  unable  to  capture  the  entire 


Fig.  4.  Model  validation  using  rat  core  body  temperature  data  obtained  from 
radiotelemetry  measurements  of  control  rats  ( n  =  30)  placed  at  an  ambient 
temperature  of  22°C  (A)  and  heat-stressed  rats  ( n  =  50)  in  an  incubator 
maintained  at  37°C  ( B )  for  163  min.  The  dashed  lines  represent  the  average  of 
the  measured  core  temperature  values  and  the  gray  bars  indicate  one  standard 
deviation.  The  core  temperatures  predicted  by  our  model  are  plotted  by  solid 
lines. 


time  evolution  of  the  experimental  core  temperature  profile. 
This  may  be  due  to  the  use  of  abrupt  lights-on  and  -off 
transitions  to  simulate  day-night  alternation  in  the  experiments, 
which  led  to  abrupt  increases  in  rat  core  temperature  in  antic¬ 
ipation  of  the  lights-off  (active)  period  at  1800.  This  response 
is  anticipatory  of  the  active  period  when  feeding,  grooming, 
and  other  behaviors  are  manifested  during  the  normal  circadian 
cycle.  The  effect  has  been  extensively  studied  in  the  past  (69, 
77),  although  in  relation  to  sleep  entrainment.  The  natural 
day-night  transition,  however,  does  not  occur  abruptly;  rather, 
the  light  intensity  changes  continuously  throughout  the  period. 
Our  computational  model  was  unable  to  mimic  the  abrupt  light 
and  dark  transitions  and  instead  incorporated  a  continuous 
day-night  transition. 

Figure  4 B  shows  the  comparison  of  predicted  and  measured 
core  temperatures  for  the  heat-stressed  animals.  Again,  the 
model  predictions  followed  the  experimental  data  ( n  =  50) 
closely,  capturing  the  magnitude  of  rise  in  core  temperature 
during  heat  stress  as  well  as  the  post-heat-stress  temperature 
response.  The  maximum  rise  in  temperature  due  to  heat  stress 
was  — 4°C.  This  rise  in  core  temperature  was  fixed  in  the 
experiments  since  animals  were  removed  from  the  incubator  as 
soon  as  they  reached  a  temperature  of  41.8°C.  The  computa¬ 
tional  model  precisely  predicted  the  timing  and  amplitude  of 
the  temperature  rise,  which  were  not  fixed  in  the  simulations 
because  the  model  input  boundary  conditions  only  included  the 
measured  time-dependent  temperature  profile  in  the  incubator, 
ranging  from  22°C  to  37°C  (per  the  experimental  protocol), 
thereby  further  validating  the  model  predictions.  After  the 
removal  of  heat  stress  (at  —1045),  the  core  temperature  started 
to  fall  rapidly  but  remained  —  1.0°C  higher  in  heat-stressed  rats 
compared  with  control  rats  even  after  110  min  of  cooling.  A 
delayed  nighttime  circadian  rise  in  temperature  was  observed 
in  the  measured  data  but  was  not  captured  by  the  simulations. 

Spatial  Temperature  Distribution 

Figure  5  shows  the  predicted  spatial  temperature  distribution 
in  control  (Fig.  5A)  and  heated  (Fig.  5 B)  rats  at  the  time  when 
the  core  temperature  reached  its  peak  value  in  the  heated 
animal.  In  Fig.  5,  A  and  B ,  different  color  scales  for  tempera¬ 
ture  were  used  to  illustrate  the  temperature  nonuniformity  in 
different  organs.  For  the  control  rat  (Fig.  5A),  the  temperature 
profile  on  the  skin  surface  (left  panel)  was  fairly  uniform 
except  in  the  tail,  where  lower  temperatures  were  observed. 
The  average  tail  surface  temperature  of  24.3°C  was  signifi¬ 
cantly  lower  than  the  average  skin  temperature  (excluding  the 
tail)  of  36.5 °C,  as  corroborated  by  previously  published  exper¬ 
imental  data  (41,  70,  78).  The  extremities,  forelimbs  at  an 
average  temperature  of  35.7°C  and  hindlimbs  at  36.0°C,  were 
also  cooler  than  other  regions  of  the  body  surface  (36.7 °C). 
The  temperature  distribution  for  the  inner  organs  (Fig.  5A,  right 
panel)  shows  comparatively  higher  temperatures  in  the  liver 
and  head  regions.  In  the  skeletal  system,  the  bones  in  the 
forelimbs  (humerus,  radius,  and  ulna)  and  those  in  the 
hindlimbs  (femur,  tibia,  and  fibula)  were  found  to  be  at  a  lower 
temperature  compared  with  the  vertebrae  and  ribs. 

For  the  heat-stressed  rat  (Fig.  5 B),  there  was  a  marked 
difference  in  temperatures  on  the  tail  surface  compared  with 
other  locations  on  the  body  surface  (left  panel),  similar  to  the 
control  rat.  However,  the  magnitude  of  temperature  difference 
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Fig.  5.  Predicted  spatial  temperature  distribution  at  time  =  4.75  h  from  the  start  of  the  simulations  (when  the  predicted  core  body  temperature  of  the  heat-stressed 
rat  reached  its  maximum  value)  in  the  control  (A)  and  heat-stressed  rats  (. B ).  In  A  and  B,  different  color  scales  were  used  to  illustrate  the  spatial  temperature 
distribution  across  the  various  organs.  Left,  skin  surface  temperatures;  right  internal  organ  temperatures. 


was  significantly  smaller,  1.7°C  as  opposed  to  12.2°C  in  the 
control  rat.  Likewise,  the  forelimbs  and  hindlimbs  were  also  at 
a  lower  temperature  than  the  other  regions  of  the  body  surface. 
The  most  critical  observation  from  the  temperature  distribution 
in  the  inner  organs  of  the  heat-stressed  rat  (Fig.  5 B,  right  panel) 
was  the  occurrence  of  considerably  higher  temperatures  in  the 
liver.  This  result  provides  an  indication  of  the  increased  sus¬ 
ceptibility  of  the  liver  to  heat  injury.  The  lower  ribs  also  had 
higher  temperatures  compared  with  other  bones  due  to  their 
proximity  to  the  liver.  Furthermore,  temperatures  in  the  head 
and  heart  were  also  higher  than  other  organs.  Although  the 
model  predictions  provide  critical  insights  about  the  expected 
temperature  distribution  in  the  different  organs  during  heat 
stress,  it  is  important  that  these  predictions  are  experimentally 
validated  in  the  future. 

Thermoeffector  Responses  to  Heat  Stress 

Figure  6  shows  the  contribution  of  the  major  heat  transfer 
mechanisms  in  the  control  and  heat-stressed  animals  at  the  time 
corresponding  to  the  observance  of  peak  core  temperature  and 
immediately  after  the  removal  from  heat  stress.  Heat  loss  due 
to  convective  heat  transfer  with  the  surroundings  decreased 
during  heat  stress  as  the  incubator  temperature  was  higher 
compared  with  the  control  ambient  condition  (37°C  vs.  22°C). 
Immediately  after  the  animal  was  placed  at  room  temperature 
(22°C),  convective  heat  loss  from  the  animal  increased  due  to 
a  larger  temperature  difference,  which  accelerated  cooling. 
Similarly,  heat  was  lost  due  to  increased  evaporation  during 
heat  stress,  which  was  sustained  even  after  heat  stress.  Heat 
generated  due  to  metabolism  increased  in  the  animal  during 
heat  stress  due  to  an  increase  in  whole  body  temperatures  (Fig. 
5 B).  Blood  flow  contributed  to  an  increase  in  heat  transfer 
during  heat  stress.  However,  after  the  removal  of  heat  stress, 
blood  circulation  facilitated  cooling  in  the  animal. 


The  Heat-Stress  Response  in  the  Organs 

Temperature  variations  in  the  organs  computed  by  the  model 
were  plotted  as  a  function  of  time  to  study  the  differential 
heat-stress  response  and  are  shown  in  Fig.  7.  The  temperatures 
shown  in  the  plot  are  volume-weighted  average  temperatures  in 
the  entire  organ  and  not  the  maximum  temperature  at  a  single 
location  inside  the  organ,  thereby  providing  an  index  for 
overall  heat  load  to  a  particular  organ.  The  highest  predicted 
temperatures  were  observed  in  the  liver,  followed  by  the  brain, 
stomach,  and  heart.  The  temperatures  in  these  organs  (solid 
lines)  were  consistently  higher  than  the  core  temperature 
(dashed  lines)  throughout  the  day-night  cycle.  The  difference 
between  the  maximum  temperatures  in  these  organs  and  the 
core  were  1.2,  0.8,  0.4,  and  0.4°C  for  the  liver,  brain,  stomach, 


8  n  □Control  HHeat  stress  ■Immediately  after  heat  stress 


-8  J 


Convection  Evaporation  Metabolism  Blood  Flow 

Fig.  6.  Contribution  of  major  heat  transfer  mechanisms  computed  by  the  model 
in  the  control  and  heat-stressed  animals  at  time  =  4.75  h  from  the  start  of  the 
simulation  (corresponding  to  the  observance  of  peak  core  temperature)  and  for 
the  time  immediately  after  heat  stress  (exposed  animals  only).  Positive  values 
indicate  heat  gained  by  the  animal,  whereas  negative  values  denote  heat  lost  by 
the  animal. 
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Fig.  7.  Differential  thermal  responses  in  various  organs  (solid  lines)  compared  with  the  core  temperature  response  (dashed  lines).  The  horizontal  dotted  line  in 
each  plot  indicates  a  temperature  of  41.8°C,  the  temperature  to  which  the  core  was  restricted  in  the  experiments.  The  plots  are  ordered  from  highest  (top)  to  lowest 
(bottom)  organ  temperature. 


and  heart,  respectively.  The  kidneys,  lungs,  and  intestines  were 
nearly  at  the  same  temperature  as  the  core,  with  the  difference 
in  peak  temperatures  being  <0.1  °C.  The  average  temperature 
in  the  tail  remained  considerably  lower  than  the  core  temper¬ 
ature  for  the  entire  duration  of  the  day-night  cycle.  Peak 
temperatures  for  all  organs  were  observed  at  the  same  time  as 
the  core  except  for  the  kidneys,  lungs,  and  intestines,  which 
were  out  of  phase  with  a  delay  of  —30  min. 

To  compare  the  thermal  response  in  heat-stressed  rats  with 
controls,  the  difference  in  average  peak  temperature  in  the 
various  organs  was  calculated  and  is  shown  in  Table  4.  Aver¬ 
age  peak  temperature  was  defined  as  the  mean  temperature  in 
an  organ  at  the  time  when  the  highest  temperature  was  ob¬ 
served  in  that  organ  during  the  heat- stress  condition.  The 
temperature  difference  in  each  organ  was  calculated  by  deter¬ 


mining  the  average  peak  temperature  in  the  heat- stress  condi¬ 
tion  and  subtracting  it  by  the  corresponding  average  tempera¬ 
ture  in  the  same  organ  at  the  same  instant  of  time  in  the  control 
rat.  The  highest  rise  in  temperature  was  observed  in  the  tail, 
which  was  more  than  three  times  the  rise  in  core  temperature 
(13.7°C  vs.  4.2°C).  There  was  also  a  significant  rise  in  tem¬ 
perature  in  the  liver  and  brain,  while  the  stomach  exhibited  an 
equal  rise  in  temperature  compared  with  the  core.  For  all  other 
organs,  the  increase  in  temperature  due  to  heat  stress  was  lower 
than  that  in  the  core.  The  lowest  increase  in  temperature  was 
observed  in  the  blood. 

To  test  the  robustness  of  our  model  predictions,  we  performed 
simulations  for  a  range  of  incubator  temperatures  (37.0  ±  1.0°C) 
in  the  vicinity  of  the  actual  experimental  temperature  (37.0°C) 
used  to  induce  heat  stress  in  the  animals.  These  simulations  also 
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Table  4.  Statistics  of  the  predicted  heat-stress  responses  in 
the  various  organs 


Location 

Average  Peak 
Temp,  °C 

Increase  in  Average 
Peak  Temp  vs. 
Control,  °C 

Maximum 
Temp,  °C 

Organ  Volume 
Fraction  with 
r>41.8°C,  % 

Liver 

42.4 

4.7 

42.9 

95.6 

Brain 

42.0 

4.6 

42.2 

91.2 

Stomach 

41.6 

4.2 

42.3 

12.1 

Heart 

41.6 

4.1 

42.2 

15.5 

Kidneys 

41.3 

3.9 

41.9 

0.2 

Core 

41.2 

4.2 

41.2 

Lungs 

41.2 

3.9 

41.6 

0.0 

Intestines 

41.1 

3.8 

42.0 

0.4 

Bone 

41.0 

3.9 

42.4 

2.6 

Blood 

40.4 

3.1 

40.4 

Tail 

38.5 

13.7 

40.7 

0.0 

Average  peak  temperature  was  defined  as  the  mean  temperature  in  an  organ 
at  the  time  when  the  highest  temperature  was  observed  in  that  organ  during 
heat  stress.  Maximum  temperature  corresponded  to  the  highest  temperature  at 
any  location  in  an  organ  during  heat  stress.  The  organs  in  the  table  are  ordered 
from  highest  (top)  to  lowest  (bottom)  organ  temperature. 

tested  any  potential  effect  that  the  precision  (±0.2°C)  of  the 
incubator  temperature  sensor  could  have  in  the  predicted  organ 
temperatures.  Figure  8  shows  organ-averaged  temperature 
changes  in  the  different  organs  relative  to  the  core  temperature  for 
incubator  temperatures  of  36.0,  36.8,  37.0,  37.2,  and  38.0°C.  The 
model  predicted  that  the  highest  temperatures  occurred  in  the  liver 
for  all  cases,  followed  by  brain  and  stomach.  In  addition,  we 
found  that  the  organs  had  the  same  heat  load  relative  to  each  other 
in  all  simulations. 

Average  temperature  rise  provides  an  index  of  the  heat- 
stress  response  for  the  overall  organ.  However,  during  heat 
stress,  there  may  be  locations  inside  the  organs  that  may  be  at 
a  considerably  higher  temperature  than  the  average  value. 
These  regions  may  be  more  susceptible  to  post-heat- stress 
disorders  compared  with  other  locations  in  the  organ  and, 
therefore,  identification  of  these  “hot  spots”  is  critical.  The 
computed  maximum  temperatures  inside  each  organ  are  shown 
in  Table  4.  The  maximum  temperatures  in  each  organ,  except 
for  the  tail,  were  higher  than  the  core  temperature.  For  exam¬ 
ple,  in  the  liver,  the  maximum  temperature  reached  2°C  higher 
than  the  measured  core  temperature.  In  addition  to  the  liver, 
brain,  heart,  and  stomach,  which  had  average  temperatures 
higher  than  the  core,  other  organs  that  had  average  tempera¬ 
tures  similar  to  that  of  the  core,  such  as  the  kidneys,  lungs, 
intestines,  and  some  bones,  also  had  hot  spots  that  exceeded  the 
observed  core  temperature.  In  certain  organs,  such  as  the 


intestine,  a  significantly  higher  maximum  temperature  com¬ 
pared  with  the  average  temperature  was  observed. 

The  portion  of  an  organ  reaching  a  temperature  greater  than 
the  lethal  temperature  for  the  species  could  also  be  an  impor¬ 
tant  quantitative  indicator  of  the  heat-stress  response.  The 
lethal  core  temperature  for  rats  has  been  reported  to  be  in  the 
range  of  42.5°C  to  44.8°C  (33).  In  this  work,  the  core  temper¬ 
ature  was  restrained  to  41.8  °C  to  limit  animal  mortality  during 
the  course  of  the  experiment.  Therefore,  to  quantify  the  heat 
load  in  each  organ,  calculations  were  performed  using  the 
computational  model  to  determine  the  volume  fraction  of  the 
organs  reaching  the  limiting  temperature  of  41.8°C  or  higher. 
The  volume  fractions  are  shown  as  percentages  in  Table  4.  As 
expected,  and  in  agreement  with  other  measures  of  the  heat- 
stress  response  discussed  above,  the  liver  and  brain  showed  the 
greatest  susceptibility  to  heat  load,  with  >90%  of  these  organs 
reaching  a  temperature  >41.8°C.  In  contrast,  10-15%  of  the 
heart  and  stomach  also  attained  temperatures  >41.8°C.  Other 
organs  had  negligible  percentages  of  their  volumes  reaching 
the  limiting  temperature. 

DISCUSSION 

In  this  study,  we  developed  an  anatomically  accurate  com¬ 
putational  model  for  thermoregulation  in  rats  during  heat  stress 
and  compared  the  heat- stress  response  with  the  control  condi¬ 
tion.  The  novelty  of  this  work  stems  from  the  comprehensive 
description  of  various  heat  transfer  mechanisms  and  incorpo¬ 
ration  of  an  accurate  representation  of  the  rat  anatomy  to 
demonstrate  that  monitoring  the  differential  thermal  response 
in  each  body  organ  is  critical  to  assess  the  actual  condition  of 
the  animal  during  heat  stress.  While  anatomically  accurate 
computational  models  have  been  developed  for  thermoregula¬ 
tion  in  rats  (44,  71,  75),  the  aims  of  these  studies  were  to 
investigate  the  effects  of  electromagnetic  radiations  over  short 
intervals  of  time  (30-5,400  s)  and,  therefore,  did  not  have  to 
consider  certain  key  phenomena  with  long  periodicity.  In 
addition,  these  studies  did  not  investigate  the  differential  ther¬ 
mal  response  in  the  organs,  which  can  potentially  lead  to 
distinct  severity  of  organ  dysfunction  during  severe  heat  stress. 
A  comprehensive  model  for  describing  the  differential  heat- 
stress  response  for  an  entire  day-night  cycle  has  not  been 
previously  reported. 

Mechanistic  models,  such  as  the  one  developed  in  this  work, 
are  models  that  incorporate  and  enable  understanding  of  the 
underlying  biological  processes  and,  therefore,  provide  predic¬ 
tive  capability  over  a  wide  range  of  conditions.  A  critical 


Fig.  8.  Organ- averaged  temperature  changes 
relative  to  the  core  temperature  for  five  sim¬ 
ulations  of  incubator  temperatures  ranging 
from  36°C  to  38°C.  Positive  values  indicate 
higher  average  organ  temperature  compared 
with  the  core,  whereas  negative  values  de¬ 
note  lower  organ  temperature. 


Liver  Brain  Stomach  Heart  Kidneys  Lungs  Intestines  Bone  Tail 
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element  in  the  development  of  mechanistic  models,  which  has 
been  lacking  in  similar  applications  (71,  75),  is  a  detailed 
model  validation  procedure  using  in  vivo  experimental  data.  In 
some  studies  (40),  the  same  set  of  experimental  data  was  used 
for  estimating  the  model  parameters  and  validating  the  model 
predictions,  which  limits  the  applicability  of  the  model  to  “fit” 
the  observed  conditions.  In  contrast,  in  our  model  validations, 
we  compared  our  model  predictions  with  extensive  indepen¬ 
dent  experimental  data,  thereby  retaining  the  mechanistic  na¬ 
ture  of  the  model  for  making  predictions  beyond  the  observed 
conditions.  The  experimental  data  used  for  model  validation 
were  not  used  to  estimate  the  model  input  parameters.  We  used 
a  large  number  of  animal  experiments  to  quantify  the  average 
thermal  response  of  both  heat-stressed  animals  ( n  =  50)  and 
control  animals  (n  =  30)  and  compared  them  to  the  model 
predictions  (Fig.  4).  Moreover,  additional  published  data  were 
used  to  corroborate  the  accuracy  of  the  model  predictions  of 
certain  physiological  variables  (Table  3). 

Physiologically,  the  temperature  of  the  blood  is  a  function  of 
the  environmental  conditions  and  heat  transfer  characteristics 
similar  to  any  other  tissue  or  organ  in  the  body.  However,  in  a 
large  number  of  computational  models  used  to  predict  temper¬ 
ature  distribution,  blood  temperature  is  assumed  to  be  a  con¬ 
stant  (23,  40,  57,  71,  75,  90).  This  assumption  may  only  be 
valid  for  cases  when  a  portion  of  the  body  is  being  heated  and 
the  overall  temperature  changes  in  the  body  may  be  considered 
insignificant.  When  the  aim  of  the  model  is  to  study  the 
thermoregulation  in  body  organs  over  a  period  of  time,  it  is 
critical  to  consider  changes  in  blood  temperature  as  a  function 
of  temperature  changes  in  other  organs.  Otherwise,  the  model 
formulation  would  artificially  introduce  a  source  (or  sink)  term 
for  generation  (or  loss)  of  heat  due  to  the  assumption  of 
constant  blood  temperature,  which  is  not  physiological.  Our 
model  accounts  for  this  change  in  blood  temperature  using  a 
novel  methodology,  which  incorporated  the  mixing  of  blood 
from  various  organs,  and  is  therefore  able  to  more  accurately 
capture  the  experimental  observations. 

Even  at  an  ambient  temperature  of  22 °C,  at  which  rats 
demonstrate  normothermic  core  temperatures  (33),  there  is 
substantial  variation  in  temperature  distribution  in  the  different 
organs  (14,  15,  34,  70).  To  accurately  study  the  resulting 
dynamics  of  evolution  of  temperature  profiles  in  the  organs  due 
to  subsequent  heat  stress,  it  is  important  to  start  with  the  correct 
normothermic  response.  Consequently,  it  is  critical  that  the 
model  simulations  of  the  steady- state  baseline  conditions  are 
capable  of  predicting  the  normothermic  temperature  distribu¬ 
tion  as  opposed  to  having  to  provide  the  starting  temperature 
distribution  as  an  initial  condition.  This  not  only  ensures  that 
the  model  is  initialized  in  a  physiologically  correct  manner, 
which  can  be  validated,  but  also  justifies  the  selection  of  input 
parameters  and  heat- transfer  mechanisms  for  the  model.  Al¬ 
though  some  of  the  published  models  have  used  steady-state 
simulations  to  initialize  the  model  and  predict  the  baseline 
response  (75,  76),  a  majority  have  used  either  a  uniform 
temperature  distribution  for  the  entire  body  as  the  initial 
condition  (40)  or  a  steady- state  simulation-based  temperature 
distribution  that  was  not  validated  (27).  In  our  modeling 
framework,  we  predicted  the  temperature  distribution  in  nor- 
mothermia  as  the  result  of  a  steady- state  simulation  and  com¬ 
pared  the  predictions  with  published  physiological  data. 


Past  studies  have  emphasized  the  need  to  characterize  heat 
transfer  in  the  rat’s  tail  in  detail  because  of  its  role  in  facili¬ 
tating  heat  loss  during  heat  stress  (67,  78,  87).  The  tail,  owing 
to  its  large  surface  area  and  an  enhanced  blood  supply,  can 
easily  exchange  heat  with  the  surroundings  by  convection  and 
radiation  (20).  The  tail,  therefore,  serves  as  an  important  organ 
in  thermoregulation,  vasodilating  blood  vessels  to  permit  in¬ 
creased  blood  perfusion  for  greater  cooling  during  heat  stress 
and  vasoconstricting  them  to  minimize  heat  loss  when  placed 
in  a  cold  environment.  Although  these  features  have  been 
implemented  in  stand-alone  models  for  tail  thermoregulation 
(67,  87),  whole  body  models  have  not  accounted  for  these 
additional  heat  transfer  mechanisms  (40,  75).  The  change  in 
heat  transfer  characteristics  in  the  tail  at  different  ambient 
temperatures  is  primarily  attributed  to  the  change  in  convective 
heat  transfer  and  rate  of  blood  perfusion.  These  mechanisms 
were  incorporated  in  our  model,  and,  consequently,  we  could 
more  accurately  capture  the  temperature  distribution  in  the  tail 
at  normal  ambient  temperature  as  well  as  during  heat  stress. 
Our  model  predicted  that  convection  was  the  dominant  mech¬ 
anism  of  heat  transfer  in  the  tail  (results  not  shown). 

In  addition  to  significant  difference  in  temperature  distribu¬ 
tion  in  the  different  organs  at  any  given  time,  the  temperatures 
in  the  organs  also  change  during  the  day-night  cycle  due  to 
circadian  rhythmicity.  Metabolic  heat  generation  and  blood 
perfusion  have  been  known  to  demonstrate  circadian  variations 
(21,  68).  This  work  is  the  first  attempt  to  predict  the  circadian 
changes  in  temperatures  using  a  mechanistic  whole  body 
model.  Our  model  accurately  captures  the  nocturnal  character¬ 
istic  in  the  circadian  variation  of  temperature  with  higher 
nighttime  temperatures  compared  with  daytime  temperatures. 
Incorporation  of  this  circadian  rhythm  is  important  in  the 
context  of  heat  stress  because,  compared  with  the  rise  in 
temperature  during  stress  (~4°C),  the  circadian  variation  in 
temperature  is  appreciable  [~0.8°C  for  this  study,  but  higher 
(1.2°C)  in  previous  reports  (33,  59,  68)].  The  extent  of  heat 
injury,  therefore,  depends  not  only  on  the  external  temperature 
but  also  on  the  time  of  day  it  is  initiated,  as  even  a  lower 
ambient  temperature  may  elicit  a  similar  thermal  response  if 
applied  at  a  different  time  during  the  day.  Preclinically,  the 
model  can  be  used  to  design  optimized  heating  experiments  to 
elicit  desired  responses  in  animals  to  study  heat-stress  progres¬ 
sion. 

The  dynamics  of  temperature  changes  during  heat  stress  is  a 
result  of  complex  interplay  of  anatomy,  tissue  properties,  and 
external  and  internal  heat- transfer  mechanisms.  As  the  ambient 
temperature  increases,  temperatures  at  the  rat’s  skin  surface 
start  to  rise  due  to  increased  convective  heating  (Fig.  4 B).  Heat 
is  subsequently  transferred  to  the  internal  organs  by  conduction 
and  convection  due  to  blood  perfusion.  Conduction  is,  in  turn, 
dependent  not  only  on  the  characteristics  of  the  particular 
tissue  but  also  on  the  anatomical  placement  of  the  organ,  and, 
therefore,  organs  farther  away  from  the  skin  surface  are  heated 
later,  generating  a  temperature  gradient  in  the  body.  Con¬ 
versely,  convective  heat  transfer  due  to  blood  flow  is  much 
faster  and  tends  to  equilibrate  the  temperatures  in  the  different 
organs.  The  rate  of  convective  heat  transfer  corresponds  to  the 
temperature-dependent  blood  perfusion  rate  in  the  particular 
organ.  Moreover,  different  organs  have  varied  metabolic  heat 
generation  rates,  which  themselves  change  due  to  change  in 
temperatures.  At  the  same  time,  heat  is  lost  from  the  skin 
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surface  due  to  evaporation  of  saliva.  All  these  mechanisms  of 
heat  transfer  make  the  process  highly  nonlinear  with  respect  to 
any  variable,  e.g.,  although  the  heart  has  the  highest  metabolic 
heat  generation  rate  (Table  1),  it  does  not  reach  the  highest 
temperature  due  to  the  combined  effect  of  all  other  heat- 
transfer  mechanisms.  The  computational  model  incorporates 
the  various  phenomena  involved  in  the  heat- stress  process 
necessary  for  predicting  the  temperature  responses,  as  shown 
in  Figs.  4,  5,  and  7,  and  enabling  us  to  tease  out  the  key 
heat- transfer  mechanisms  (Fig.  6).  The  temperature  rise  in 
most  organs  is  in  phase  (Fig.  7),  possibly  because  heat  transfer 
due  to  blood  perfusion  is  the  dominant  mode  of  heating.  Heat 
conduction  also  plays  an  important  role,  as  can  be  observed  by 
the  temperature  variations  within  an  organ  (Fig.  5 B  and  Table 
4).  For  organs,  such  as  lungs  and  kidneys,  which  reach  the 
maximum  temperature  after  a  delay  (Fig.  7),  conduction  is  the 
more  significant  heat-transfer  mechanism. 

The  most  powerful  contribution  of  the  model  is  its  capability 
to  provide  a  framework  for  comprehensive  whole  body  pre¬ 
diction  of  temperature  changes.  The  model  results  show  that 
certain  organs  attain  much  higher  peak  temperatures  compared 
with  the  core  (Table  4).  Measuring  only  the  core  or  rectal 
temperature,  as  routinely  done  in  heat- stress  studies  (7,  54,  88), 
may  not  present  an  accurate  picture  of  the  thermal  response  in 
the  entire  body.  Additionally,  significant  differences  in  tem¬ 
perature  distribution  can  be  observed  even  for  organs  in  close 
proximity  to  each  other  (Figure  5 B  and  Table  4).  Therefore, 
irrespective  of  where  the  temperature  measurement  is  done,  it 
is  critical  to  consider  the  spatial  variation  of  temperature  in 
other  organs.  Another  drawback  of  measuring  only  the  core 
temperature,  or  the  temperature  at  a  specific  location  in  an 
organ,  is  that  it  does  not  provide  insights  about  the  possible 
presence  of  “hot  spots”  or  regions  of  high  temperatures  within 
an  organ.  Our  model  demonstrates  (Table  4)  that  there  may  be 
locations  within  an  organ  that  may  achieve  considerably  high 
temperatures  even  though  other  regions  in  the  organ  may  be 
well  below  the  lethal  temperature.  This,  in  turn,  makes  analysis 
of  spatial  temperature  distribution  all  the  more  essential.  Al¬ 
though  imaging  techniques,  such  as  MRI-based  thermometry, 
have  been  used  for  two-dimensional  spatial  monitoring  of 
temperature  changes  in  specific  organs  during  thermal  therapy 
of  short  durations  (2,  66),  live  monitoring  of  3-D  temperature 
distributions  in  the  entire  body  in  large  groups  of  animals  in  a 
heat-stress  study  is  still  not  feasible.  Therefore,  with  the  limited 
availability  of  experimental  methods,  our  computational  model 
nicely  complements  core  temperature  measurement  experi¬ 
ments  in  providing  an  estimate  of  the  temperature  distribution 
in  the  entire  body  during  heat  stress. 

The  most  critical  finding  from  our  work  is  that,  among  the 
organs,  the  liver  experiences  the  highest  temperatures  during 
heat  stress  (Table  4).  The  brain  reached  the  second  highest 
temperatures  following  the  liver.  Based  on  the  expected  tem¬ 
perature  elevation  in  the  various  organs  predicted  by  our 
model,  appropriate  cooling  strategies  could  be  used  for  specific 
organs.  Furthermore,  histopathology  studies  of  heat-stressed 
animals  have  shown  progressively  greater  organ  damage  asso¬ 
ciated  with  higher  temperatures  (56).  Liver  dysfunction  has 
indeed  been  reported  to  occur  early  in  heat  stroke  (31,  83). 
Moreover,  cerebral  ischemia,  injury,  and  dysfunction  during 
heat  stroke  have  also  been  reported  in  the  literature  (12,  13). 
While  further  experiments  and  analyses  are  needed  to  directly 


correlate  the  subsequent  inflammatory  response  and  MODS 
with  temperature  distribution  in  the  organs,  the  model  provides 
a  reasonable  estimate  of  the  susceptibility  of  a  particular  organ 
for  dysfunction  as  a  result  of  heat  stress. 

Although  our  model  includes  a  comprehensive  array  of  heat- 
transfer  mechanisms  involved  during  thermal  stress  in  rat,  ther¬ 
moregulation  is  a  complicated  process,  and  it  is  not  possible  to 
replicate  all  mechanisms  involved  in  a  living  organism.  Another 
limitation  is  that  assumptions  were  made  throughout  the  model 
formulation  stage.  For  example,  heat  transfer  due  to  blood  perfu¬ 
sion  was  simplified  by  using  the  Pennes  bioheat  transfer  equation 
and  by  employing  a  mixing  formulation  for  change  in  blood 
temperature.  Nevertheless,  the  model  was  extensively  validated 
using  in  vivo  experimental  data  to  justify  these  and  other  model 
assumptions  as  well  as  the  selection  of  the  key  heat-transfer 
mechanisms.  Moreover,  no  ad  hoc  model  input  parameters  were 
used,  so  the  model  remains  mechanistic  in  nature. 

The  progression  of  heat  stroke  from  the  temperature  rise  in 
the  body  to  inflammation  and  MODS  is  a  complex  process. 
Our  model  aims  at  investigating  only  the  initial  phase  of  this 
cascade  of  events.  However,  to  model  the  entire  pathophysi¬ 
ology,  it  is  important  not  only  to  include  the  physiological 
events  (temperature  change)  but  also  the  resulting  biological 
response  (inflammation).  For  example,  the  feverlike  response 
in  core  temperature  after  an  initial  temperature  rise  cannot  be 
predicted  mechanistically  by  the  model  without  incorporating 
the  effects  of  biological  factors,  such  as  the  role  of  cytokines 
(53),  activated  as  a  result  of  heat  injury.  However,  as  with  any 
computational  model,  to  incorporate  these  relationships  in  a 
model,  it  is  important  to  first  study  them  experimentally.  A 
detailed  experimental  animal  model  is,  therefore,  needed  to 
clearly  elucidate  the  relationship  between  the  temperature  rise 
in  each  organ  and  the  resulting  change  in  inflammatory  medi¬ 
ators  before  such  a  molecular  computational  model  can  be 
envisioned.  Nevertheless,  our  present  model  provides  a  starting 
point  for  the  development  of  such  a  multiscale,  molecular- 
physiological  computational  model  for  heat  stroke. 

Monitoring  of  organ  temperatures  is  not  a  trivial  endeavor  in 
a  conscious  rodent  model  because  multiple  transmitters  cannot 
be  implanted  into  the  peritoneal  cavity  of  an  animal  (large  size 
precludes  this)  and  regional  differences  within  an  organ  would 
need  to  be  considered  as  well.  However,  brain  (e.g.,  hypotha¬ 
lamic,  cerebellum)  temperature  could  be  simultaneously  mon¬ 
itored  with  the  core  temperature.  It  might  also  be  possible  to 
monitor  individual  organ  temperatures  under  anesthetized  con¬ 
ditions.  Presumably,  with  advancements  in  radiotelemetry  (i.e., 
smaller  transmitter  size  and  elimination  of  frequency  interfer¬ 
ence  among  devices)  and  imaging  technologies,  it  may  be 
possible  to  validate  organ- specific  temperature  predictions  in 
the  near  future.  Alternatively,  one  could  take  a  more  simple 
approach  of  rapidly  sacrificing  animals  and  measuring  organ 
temperature  in  situ  at  different  times  during  recovery  to  vali¬ 
date  the  model  predictions,  although  with  limited  accuracy. 

The  ability  to  predict  the  spatio-temporal  temperature  dis¬ 
tribution  opens  opportunities  for  potential  clinical  and  design 
applications.  Because  the  model  is  mechanistic,  it  can  be  used 
for  simulating  a  multitude  of  “what-if  ’  scenarios  of  different 
heat-stress  episodes  to  design  optimized  experimental  proto¬ 
cols.  This  would  help  in  substantially  reducing  the  number  of 
animals  required  for  such  experimentations.  The  model  could 
also  be  used  to  test  and  design  novel  heat  stress  management 
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techniques,  such  as  localized  cooling.  Brain  cooling  has  shown 
promise  in  the  recovery  from  heat  stroke  in  rats  (47),  and 
similar  localized  cooling  techniques  could  be  used  for  inducing 
hypothermia  in  other  heat- stroke- susceptible  organs,  such  as 
the  liver  (31,  83).  Further,  the  spatial  temperature  distribution 
in  the  various  organs  could  be  used  to  anchor  and  analyze 
molecular  mediators  of  heat  injury  in  an  organ- specific  man¬ 
ner,  in  particular  in  relation  to  hot  spots.  This,  in  turn,  would 
help  in  the  identification  of  predictive  and  diagnostic  biomark¬ 
ers  of  systemic  inflammation  and  organ  injury  after  heat  stress. 
If  imaging  techniques,  such  as  MRI,  are  used  to  estimate  tissue 
damage  during  heat-stress  recovery,  it  may  be  possible  to 
determine  the  effect  of  different  predicted  organ  temperatures  on 
the  progression  of  organ  damage.  Finally,  the  model  could  also  be 
used  to  simulate  the  efficacy  of  preventative  measures,  such  as 
customized  apparel,  for  limiting  heat-stress  occurrences. 

In  summary,  our  model  is  successfully  able  to  predict  the 
thermal  response  in  rats  during  heat  stress.  We  used  the  model 
results  to  demonstrate  that  monitoring  the  spatial  distribution 
of  temperature  during  heat  stress  is  critical.  The  model  results 
can  be  used  to  quantify  the  heat  load  in  various  organs  by 
determining  indices,  such  as  average  temperature  rise,  maxi¬ 
mum  temperature,  and  volume  fraction  of  organ  reaching  the 
fatal  temperature,  which  may  not  be  possible  to  experimentally 
measure.  The  model,  therefore,  complements  experiments  and 
can  be  used  for  further  understanding  the  complex  processes 
involved  in  thermoregulation  in  rats.  Finally,  the  model  enables 
identification  of  organ- specific  risks  during  heat  stress,  which 
will  potentially  aid  in  the  development  of  improved  clinical 
strategies  for  injury  prevention  and  management. 
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